Nonexponential decoherence and momentum subdiffusion 
in a quantum Levy kicked rotator 



O 
O 

> 
O 

in 



> 

(N 



o ■ 

^ : 

9 L,: 

-I— > ' 

G ■ 
cd 

3 : 



Henning Schomerus 1 and Eric Lutz 2 

1 Physics Department, Lancaster University, Lancaster, LAI 4YB, UK 
2 Department of Physics, University of Augsburg, D-86135 Augsburg, Germany 

(Dated: February 1, 2008) 

We investigate decoherence in the quantum kicked rotator (modelling cold atoms in a pulsed 
optical field) subjected to noise with power-law tail waiting-time distributions of variable exponent 
(Levy noise). We demonstrate the existence of a regime of nonexponential decoherence where the 
notion of a decoherence rate is ill-defined. In this regime, dynamical localization is never fully 
destroyed, indicating that the dynamics of the quantum system never reaches the classical limit. 
We show that this leads to quantum subdiffusion of the momentum, which should be observable in 
an experiment. 
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The transition between quantum and classical dynam- 
ics is not yet fully understood and the boundary be- 
tween the two remains fuzzy. Deep insight into the emer- 
gence of classical properties from quantum mechanics is 
provided by the theory of decoherence [J 0. In this 
framework, the unavoidable coupling to an external en- 
vironment induces a dynamical suppression of interfer- 
ence effects between states of a quantum system, and 
thus classical behavior. The disappearance of quantum 
interference can often be described by a decoherence fac- 
tor of the form T>(t) = exp[— inn (Ax) 2 t], where m is 
the mass of the particle, k the coupling strength and 
Ax the spatial separation of the interfering states. The 
functional form of the decoherence factor, in particu- 
lar its exponential time dependence, has been verified 
in many experiments 0, 0j @> 0] • The boundary between 
the quantum and the classical world can accordingly be 
explored by properly tuning the value of the coherence 
time, t c = l/[mn (Ax) 2 ] . This has been performed in var- 
ious ways, either by changing the separation Ax as in the 
QED cavity and ion trap experiments by the groups of 
Haroche [|| and Wineland [1] , or by modifying the value 
of the coupling constant k as done in the atom-optics 
experiments by the groups of Christensen [f| and Raizen 
Q . In yet another approach, the quantum-classical bor- 
der has been investigated by increasing the mass m of the 
system as in the double-slit experiments by the group of 
Zeilinger Q. 

Common to all these studies is the exponential func- 
tional dependence of the decoherence factor. Our aim in 
this paper is to widen our understanding of the quantum- 
to-classical crossover by proposing an environmental cou- 
pling scheme for the atom-optical experiments [B|, Q 
which allows the controlled tuning of the functional form 
of T>(t) itself, in particular of its time dependence, and 
not just of its parameters. In this way, we are able to 
uncover some hidden aspects of decoherence, especially 
of its interplay with complex quantum dynamics Q . 

A paradigm of complex quantum dynamics are cold 



atoms exposed to pulsed optical fields 0]. For chaotic 
conditions the classical motion in these systems ex- 
hibits momentum diffusion which however is suppressed 
in th e q uantum regime because of dynamical localiza- 
tion [lCj . In the recent decoherence experiments 0, [f| , 
the quantum-classical transition has been investigated 
by varying the amplitude of noisy pulses acting on the 
atoms [coupling strength k in T>(t)]. This then results 
in a non-vanishing quantum momentum diffusion with 
a renormalizcd diffusion constant. By contrast, we here 
propose to modify the length t of the time intervals be- 
tween the noisy pulses Specifically, we shall con- 
sider that the intervals are generated in a renewal process 
with a waiting-time distribution w(t) that asymptoti- 
cally behaves as a power law, w(t) oc t -1 ~ q . Power-law 
tail waiting-time distributions of this kind - generally 
referred to as Levy statistics [l2[ - occur naturally in 
many physical systems: Two prominent examples being 
subrecoil laser cooling [l3| and fluorescent intermittency 
of nanocrystal quantum dots [3| (in both these systems, 
a is close to one half). In the atom-optical setting, the 
statistics arc under the control of the experimentalist. By 
tuning the value of the exponent a, an additional param- 
eter not present in former decoherence studies, one is able 
to control the duration between successive random kicks 
given by the mean waiting time, r = / °° oIttw{t). This 
offers a unique opportunity to smoothly interpolate be- 
tween a fully coupled situation with many kicks (large a) 
to an almost isolated situation with very few kicks (small 
a) . Levy kicked systems have lately been investigated at 
the classical [l5[ as well as the quantum-mechanical level 
flC^ and unusual properties such as nonexponential be- 
havior and aging have been found. 

Interestingly, for a < 1, r becomes mathematically 
infinite. We show that this divergence of the mean wait- 
ing time dramatically affects the loss of phase-coherence 
of the atoms: the decoherence factor stops being expo- 
nential in time (instead it is given by a Mittag-Lcfflcr 
function) and the concept of a decoherence rate becomes 
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ill-defined. This results in quantum sub-diffusion of the 
momentum of the Levy kicked rotator, a unique signature 
which would be directly accessible in experiment. 

Model. The motion of an atom in a pulsed optical field 
can be modelled as a kicked rotator, a connection that 
has been established extensively in both theoretical and 
experimental works [^, [13, EH • The kicked rotator can 
be simply seen as a particle moving on a ring and period- 
ically kicked in time. In conveniently scaled coordinates, 
the Hamiltonian is given by 

2 00 

H =Y + K ncos66(t-n + + ), (1) 

n— — oo 

where the kicking potential depends on the 27r-periodic 
rotation angle 6. We assume that the particle is ini- 
tially prepared in the zero-momentum state ip(0) = \p = 
0). The subsequent stroboscopic dynamics ip(t + 1) = 
F(K t )ip(t) are generated by the Floquet operator 

F{K t ) = exp{-ih~ 1 p 2 /2)exp(~ih- 1 K t coae). (2) 

This operator describes a kick in which momentum 
changes by K t sin#, followed by a free rotation in which 
the rotation angle 9 increases by p. 

For constant (noiseless) kicks of strength K n = K > 5, 
the classical dynamics is chaotic, and on average dif- 
fusive in momentum direction: varp(i) ~ D c \t, where 
D c \ ~ K 2 /2. In the quantum regime, the kicked rota- 
tor exhibits dynamical localization [1 01 ] . an interference 
phenomenon which manifests itself in an exponentially 
decaying envelope of the quasienergy eigenstates of F(k) 
in the momentum representation. As a result of dynam- 
ical localization, momentum diffusion is suppressed and 
the variance 

varpo(*) ^ D*t*[l -exp(-t/t*)\ (with D* ~ £) cl ) (3) 

saturates after the quantum break time t* ~ D* / h 2 [l8| . 

The dramatic difference between classical and quan- 
tum momentum dynamics makes the kicked rotator ide- 
ally suited for the investigation of the crossover between 
both regimes. The crossover can be induced by sub- 
jecting the particle to additional random kicks - sim- 
ulating in such a way the noisy coupling to an exter- 
nal environment 0, Q. In a noisy kick, the kicking 
strength is slightly perturbed away from the mean value 
K, K n = K + k n , where the perturbations k n are random 
numbers with average k n = and variance k\ — k. The 
parameter k defines the strength of the noise. The time 
intervals between subsequent kicks are generated by a re- 
newal process with waiting distribution u>(t). Unlike the 
familiar Gaussian white noise, Levy noise with a power- 
law waiting time distribution is a nonstationary process 
with an autocorrelation function that depends explicitly 
on the two time variables. 

The effect of noise on the quantum kicked rotator is 
best understood using the formalism developed by Ott, 



Antonsen and Hanson In the absence of noise the 

quantum system is fully coherent: the quasienergy eigen- 
states of the Floquet operator F(K ) are exponentially lo- 
calized in momentum space, chaotic diffusion is strongly 
suppressed and the momentum diffusion constant asymp- 
totically vanishes. The effect of the external noise is to 
couple the quasienergy states and to induce transitions 
between them. As a result, dynamical localization is ex- 
tenuated and quantum diffusion takes place. There is 
thus an intimate connection between quantum momen- 
tum diffusion and phase coherence in the kicked rotator. 

A systematic, perturbative approach for calculating 
the quantum diffusion constant for weak stationary noise 
with arbitrary correlations has been put forward by Co- 
hen [13]. In this case, the survival probability of the 
quasienergy eigenstates is found to decay exponentially 
over time, D(t) = exp(—t/t c ). In the limit of weak noise, 
k <C K^/h 4 , the coherence time is given by t c ~ 2H 2 / k. 
The resulting momentum spreading is accordingly [201 ] 

varp w — ° . t + - — D -. — [1 - exp(-t/t* - t/t c )]. 
1 + ic/i* (l + <*Ac) 2 

(4) 

Hence, for stationary noise when quantum coherence is 
lost exponentially, the momentum diffuses asymptoti- 
cally with a renormalized diffusion constant D*/(l + 
t c /t*). Equation ([4]) has been successfully applied to 
quantify the decoherence process in the experiment [6|]. 

Results. In the following, we extend the calculation of 
the quantum diffusion constant beyond the perturbative 
regime: for the present case of uncorrelated perturba- 
tions K n , a random-phase approximation allows to study 
decoherence induced by strong nonstationary noise (de- 
tails of the derivation will be presented at the end of 
this paper). In a general renewal process, the number of 
noisy events N(t',t") within an interval [t',t"] (techni- 
cally known as the inverse time of the process) is a sta- 
tistical quantity. Within the random-phase approxima- 
tion, the survival probability in a quasienergy eigenstate 
over this time interval is simply given by the decoherence 
factor 

V(t', t") = exp[-N(t',t")/t c ], (5) 

where t c is the coherence time of the corresponding sta- 
tionary noise process with the same noise intensity. The 
variance of the momentum can then be obtained from 
the following expression (see details below): 

varp(i) = varpo(t) V(t, 0) + ^N{t,0) 

+ / ds varp (i — s) d s V(t, s) 
Jo 

— / ds varpo(s) d s T>(s, 0) 
Jo 

-f ds' f ds"va,rp (s' -s")d s ,d s „V(s',s").(6) 
Jo Jo 
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Here varpo(i) is the variance of the momentum in absence 
of the noise. The second term on the right-hand side 
becomes important for large noise strength. 

Stationary noise can be regarded as a special case of 
a renewal process with w{t) = 8\ tT . In the latter situ- 
ation, the number of noisy events is simply N(t',t") — 
\t' — t"\ and the decoherence function decays exponen- 
tially T>(t',t") = exp(-|t'-£"|/t c ). Equation @ then 
provides a generalization of the perturbative result (j4|) 
which accounts exactly for all correlations of the noiseless 
dynamics. Expression Q is recovered when the noiseless 
dynamics is approximated by the variance ([3|). 

We now evaluate the decoherence factor ([5]) for a re- 
newal process with a general waiting time distribution. 
The Laplace transform of V(t, 0) with respect to t reads 



a=2 1.2 1.0 0.9 



V(u,0) = 



1 



U(l + /(u)/i c )' 



(7) 



where /(it) = w(u)/(l — w(u)) is the Laplace transform 
of the sprinkling distribution f(t) = d t N(t,t') [f(t) is 
simply the probability that there is a noise event at time 
t]. For two time arguments, we have (see also Ref. [l6[ 
for similar results), 



V(t',t") =v(t',o)+t- 



dsf(s)V(t'-s,0). (8) 



Two cases have to be distinguished: For Levy noise 
with an exponent a > 1, the sprinkling distribution ap- 
proaches a constant for long times, f(t) ~ t _1 , and the 
decoherence function V(t', t") ~ exp[— \t' — t"\/(rt c )} has 
still an asymptotic exponential form. We thus have the 
important result that the effective coherence time ft c is 
directly proportional to the mean waiting time of the 
noisy kicks and therefore increases with decreasing ex- 
ponent a. In particular, the effective coherence time 
becomes infinitely large when a drops below unity, in- 
dicating that the classical limit cannot be attained in 
finite times. Clearly, the functional form of the decoher- 
ence factor has to change as well. In order to determine 
T>(t',t") for Levy noise with a < 1, we use the mathe- 
matical equivalence of Eq. {5| to the moment-generating 
function M(£) = exp(£iV) of the inverse time N, evalu- 
ated at f = -l/t c UM- We then find 



V(t,0) ^E a (-T(l + a)N(t,0)/t ( 



(9) 



l /T(an + 1) is the Mittag-Lefner 



where E a (z) 

function [22j. Asymptotically, N(t,0) ~ t a sin(7ra)/(7rc). 
The initial decay of the decoherence function is therefore 
a stretched exponential, T>(t, 0) ~ exp{i Q /[r(— a)ct c ]}, 
while for large times it crosses over to the power law 
T>(t, 0) ~ (ct c /a)t~ a . The functional dependence of the 
decoherence factor can thus be changed in a controlled 
manner by tuning the value of the exponent a. The con- 
sequences of such a modification can be directly detected 




FIG. 1: (color online). Time dependence of the momentum 
spreading varp(t) for kicked rotators which are subjected to 
Levy noise generated by a Yule-Simon distribution. The re- 
sults of numerical simulations (solid curves) are compared 
to the theoretical predictions from Eq. (|6} (dashed curves). 
In the upper panel, the exponent a varies while the noise 
strength « = 1/300 is fixed. The bottom panel shows results 
for fixed a — 0.5 and varying noise strength k. 



in the momentum diffusion of the kicked rotator: As a 
approaches unity from above, the decoherence time in- 
creases, and the asymptotic quantum diffusion coefficient 
D*/(l + rt c /t*) vanishes at a = 1. For a < 1 it follows 
from Eq. © that the momentum spreads subdiffusively, 



varp(i) 



varpo(oo) 



tr 



N(t,0) 



D*t* sin7ra 



tr. 



-t" 



(10) 



The full time dependence of the momentum spreading 
in kicked rotators with specific implementations of Levy 
noise is shown in Fig. [TJ The regular kicking strength 
is set to K = 7.5, while h = 2ir x 577/13872. The ran- 
dom numbers k n are taken from a uniform box distri- 
bution over an interval (—W,W), so that k = W 2 /3. 
The waiting-time distribution is of the Yule-Simon form 
uo(t) = «r(T)r(ct + l)/r(r + a + 1), with asymptotic 
behavior oj(t) ~ aT(a + l)/r a+1 and a mean waiting 
time t — a/ (a — 1) 23]. In the figure, the solid lines are 
computed by extensive numerical simulations of the dy- 
namics of kicked rotators over 10 5 time steps, averaged 
over 100-1000 realizations of the noise. The dashed lines 
are evaluated from Eq. (|6|), where varpo(t) is approxi- 
mated by Eq. © and D* = h 2 t* = 45.28 is obtained 
by a single-parameter fit to the results of the numeri- 
cal simulations without noise. The decoherence function 
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([5]) for the Yule-Simon distribution is determined setting 
t c = 2h 2 / k. We find good agreement between theory and 
numerics, both as a function of the exponent a as well as 
a function of the noise strength n. In particular, the nu- 
merics confirm the subdiffusive long-time dynamics (|10[) 
for a < 1, as well as the diffusive behavior for a > 1. 

Details of the calculation. The derivation of our main 
result ((6]) along with Eq. ([5|) is based on the observation 
that the phase of the transition amplitudes {r\F(Ki)\s) = 
F r ^ s (Ki) between quasienergy eigenstates \r) and \s) de- 
pends strongly on the random detuning ki = Ki — K of 
the kick. In the presence of noise we set 



F r , s {Ki) = 5 r , s exp[zuv - l/(2t c )] 



(11) 



while in the absence of noise F r ^ s (K) = <5,- iS exp(iu> r ), 
where uj r is the quasienergy of the state. 

The transition amplitudes F r<s enter the variance of the 
momentum varp(i) = t „_ C(t' ,t") when the force- 
force correlation C(t',t") — (K t iK t n sin# t / sin 9*» ) is ex- 
panded in the basis of quasienergy eigenstates [2fJ, 

C(t',i") = {r v \wa.6\s v )(a V '\wxB\r V ') x 

i'-l 

x II K «Kt»F; i+liri {Ki)F Sm+uam (K m ). (12) 

Lm—t" 



In the sum (jT3J) , only those terms with satisfy either (a) 
r n = r and s n = s, or (b) r n = s n survive the average in 
the random-phase approximation pip . The terms (b), 
however, have a zero contribution since they are multi- 
plied by (r|sin0|r) = 0, a property which follows from 
the symmetry (8,p) — > (—9,—p) of the kicked rotator. 
Each of the terms (a) further carries a weight factor 
exp[— N(t',t")/t c ], which has to be averaged using Eq. 
(|11[) . The averaging procedure over the renewal process 
then leads to the decoherence factor T>(t',t") defined in 
Eq. ([5]). Moreover, for t' = t" the noise also gives rise to 
a contribution to the typical force acting on the particle 
of the form ^"(sinfl 2 ) = K 2 /2 + («/2)/(t')< Collecting 
all terms, this yields 



C(t',t") = C (t',t")V(t',t") + -f(t')S tlA „ 



(13) 



where Co(t',t") is the force-force correlation function in 
the absence of noise. Equation ((6]) then follows by con- 
verting sums into integrals and integrating by parts. 

Conclusions. Noncxponential relaxation is familiar 
from the physics of complex systems - such as disordered 
crystals and glassy materials - and occurs whenever there 
is no clear se par ation of time scales between system and 
environment [241 ]. At the quantum level, the absence of a 
scale separation manifests itself in a noncxponential loss 
of phase coherence. We have shown that such a behavior 
can be induced in atom-optical systems by properly engi- 
neering the environment. By applying Levy noise with a 



controllable exponent a, the environment time scale - the 
mean waiting time r between noisy events - is finite for 
a > 1 and can be made divergent by choosing a < 1 . For 
the quantum Levy kicked rotator, the power-law decay 
of the decoherence factor and the nonexistence of a finite 
decoherence rate reveals itself in quantum subdiffusion 
of the momentum, a distinct signature that can easily 
be observed experimentally using cold atoms in pulsed 
optical fields. The high degree of control and tunability 
of atom-optical systems make the Levy kicked rotator an 
ideal tool for the study of the crossover from fast to slow 
decoherence. Complex systems with algebraic relaxation 
are known to never attain equilibrium; the possibility of 
having noisy quantum systems that never become classi- 
cal appears as a fascinating prospect. 
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